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We perform global fits to the parameters of the Constrained Minimal Supersymmetric Standard 
Model (CMSSM) and to a variant with non-universal Higgs masses (NUHMl). In addition to 
constraints from low-energy precision observables and the cosmological dark matter density, we 
take into account the LHC exclusions from searches in jets plus missing transverse energy signa- 
tures with about 5 fb^' of integrated luminosity at ^/s ~ 7 TeV. We also include the most recent 
upper bound on the branching ratio Bs jJ-jJ. from LHCb. The best fit of the CMSSM prefers a 
hght Higgs boson just above the experimentally excluded mass. We find that the description of 
the low-energy observables, {g — 2)^ in particular, and the non-observation of SUSY at the LHC 
become more and more incompatible within the CMSSM. A potential SM-like Higgs boson with 
mass around 126 GeV can barely be accommodated. Values for BF{Bs — > /i/i) just around the 
Standard Model prediction are naturally expected in the best fit region. 
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1. Introduction 

Current data being insufficient to constraint a general supersymmetric model, the most widely 
considered constrained model is the constrained MSSM (CMSSM) with only 5 new free parameters 
beyond the SM: Mo,Mi/2,^o,tan/3,5g«(/i), denoting respectively the universal soft supersymme- 
try breaking scalar and gaugino masses at the unification scale, the universal soft supersymmetry 
breaking trilinear scalar coupling, the ratio of the vacuum expectation values of the two CP-even 
neutral Higgs fields and the Higgs mixing parameter in the superpotential. In a minimal non- 
universal Higgs mass model (NUHMl) a universal scalar Higgs mass parameter Mh at the 
unification scale, Mhu = M^d = is added. In this paper, we investigate the global interpreta- 
tion of all existing data using our framework Fittino 



2. The Fittino framework 

The Fittino |^ ^] framework is used to perform a global Markov Chain Monte Carlo (MCMC) 
scan of the supersymmetric parameter space in all parameter dimensions. In the frequentist inter- 
pretation of the MCMC fit, assuming a Gaussian likelihood the two-dimensional la (2a) bound- 
aries are defined by A;^^ < 2.33(5.99), where A;^^ is calculated for each point by regard to the best 
fit point with the smallest x^- Using a self-optimizing chain, at least 3 million points are obtained 
within < 5.99 from the minimum for each individual fit. The SUSY particle spectrum is cal- 
culated with SPheno 3.1.4 [^, and then used in micrOMEGAs 2.2 [Q] for the prediction of the dark 
matter reUc density, in FeynHiggs 2.8.2 for the prediction of the electroweak observables, for 
the Higgs masses and for the anomalous magnetic moment of the muon , in SuperlSO 3.1 [0] 
for the flavor physics observables, and in AstroFit for the evaluation of the direct and indirect 
detection of dark matter observables. 



3. Experimental constraints 

A detained description of the experimental constraints can be found in [^]. The present and 
potential experimental measurements used in this study are indirect constraints through supersym- 
metric loop corrections and constraints from astrophysical observations. The available limits on 
SM and non-SM Higgs bosons, including the ones presented by the LHC and Tevatron collabo- 



rations at the Spring Conferences 2011, are evaluated using HiggsBounds 3.2 [10]. An extensive 
database of relevant astrophysical data are added by Astrofit to the Fittino fit process. The most 
conservative chargino mass limit of m^± > 102.5 GeV, from the direct search at LEP was included, 
leading indirectly to the exclusion of light neutralinos m o < 50 GeV in the constrained models 



considered here [|11[]. The most stringent LHC limits from the direct searches in channels with jets 
and missing transverse energy are included by emulating the search analysis. The ATLAS analy- 
sis [O] is reproduced by simulating the production of gluino and squarks with SPheno 3.1.0 and 



Herwig 2.4.2 [|13||, together with the fast detector simulation DELPHES 1.9 []14p. It is indeed not 
sufficient to only consider the 95%CL bounds published by the experimental collaborations for spe- 
cific models and particular choices of parameters. The published ATLAS limits at L,„f = 165/ pb 
and Lint = 1//^ could be precisely reproduced (see Fig. |T]), the signal grid was adapted to the AT- 
LAS analysis for L,„, = A.l / fb by reducing the systematic uncertainties. The impact of fixing Aq 
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Figure 1: Left: x contribution from the LHC SUSY search implementation compared to the pubhshed 
ATLAS and CMS Hmits ^ Right: Simulated signal yields for a different point in the (Mo,Mi/2) 

parameter space, with the SM background (gray), the CMSSM signal expectation for our grid (yellow) and 
the systematic uncertainty of 30% (orange). 

and tan j8 on the grid was checked by simulating the signal for various values of Aq and tan j3 in 
different regions of the (Mo,Mj/2) parameter space, and was found to be negligible (see Fig. |l|). 

4. Results 

The results of our CMSSM and NUHMl fits for various sets of input observables are given 
in the Table [l], for all fits we require the lightest neutralino to be the LSP, consistent radiative 
electroweak symmetry breaking and the absence of tachyons. In the plots of fitted parameters, 
the best-fit point is mai^ked by a star; all hidden dimensions having been profiled. Excluding the 
direct searches of super-symmetry at the LHC, the fitted parameter space points to a light sparticle 
spectrum, below 1 TeV (see Fig. ^)). The focus point, i.e. Mq ~ 2 TeV and M\2 ~ 150 GeV, 
is allowed in the 2a region, due to the high Higgs mass rtif, prefered by cunent data. Including 
the direct search at the LHC decreases the goodness of the fit, which arises from the coupling of 
colored and non colored sectors in the constrained models. The low energy observables are indeed 
mainly driven by non colored sparticles, whilst the channels used for the direct search at the LHC 
rely mostly on the colored one. This results in a shift upwards of tan j3 and of the masses of squarks 
and gluinos (see Fig. ||), the shift in tan/3 being cause by the correlation with the masses thi^ough 
the muon anomaly: an increase in the mass being compensated by a larger coupling in order to 
match the experimental high value of a^. Contrary to the CMSSM, the NUHMl can accomodate 
a Higgs boson as heavy as 126 GeV, depending on the value of BF{Bs — pL^pL^) (see Fig. 
The fit of the NUHMl model, including all experimental constraints, results in a large 2a contour 
with a prefered region at lower mass and the focus point to be excluded (see Fig. |^. Despite 
a lower fit tension than for the MSSM, a tension remains due to the strong coixelation between 
BF{Bs ^ H+H^), and m/, in the NUHMl. 

When the constraint of a Higgs mass of m/, = 126 GeV is included, the prefered values of 
masses and tan/3 increases (see Fig. ^), while the quality of the fit decreases further, leading to a 
tension in the fit of constrained models. This tension can hardly be relieved by leaving the top mass 
free in the fit. 
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Figure 2: Left: Parameter distributions for the LEO fit of the CMSSM with the 1 -dimensional 1(7 in red 
and the 2-dimensional 2a in blue, and the best fit point marked by a star. Right: predicted distribution of 
sparticle and Higgs boson masses from the LEO fit of the CMSSM. The full uncertainty band gives the 
1-dimensional 2a uncertainty of each mass Ax^ < 4. 
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Figure 3: Left: Parameter distributions for the LHC fit of the CMSSM with the 1-dimensional la in red 
and the 2-dimensional 2a in blue, and the best fit point marked by a star Right: predicted distribution of 
sparticle and Higgs boson masses from the LHC fit of the CMSSM. The full uncertainty band gives the 
1-dimensional 2a uncertainty of each mass Aj^ < 4. 




Figure 4: Left: Parameter distributions for the LHC+m/, fit of the CMSSM with the 1-dimensional la in 
red and the 2-dimensional 2a in blue, and the best fit point marked by a star Right: The dependence of the 
minimal "x} of the fit on m/, for different input observable sets and for the CMSSM and NUHMl. 
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Table 1: Summary of the results for the CMSSM and NUHMl fits with different sets of input observables: 
"LEO" refers to all low energy observables, "LHC" includes also the direct search for sparticles at the LHC, 
and "LHC+OT/," adds the constraint of a Higgs boson of 126 GeV. 
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Figure 5: Predicted 2a ranges of Higgs branching fractions and ratios for the LHC fit of the CMSSM with 
nil, = (126 ± 2 ± 3) GeV, including the theoretical uncertainty of the Higgs mass of 3 GeV. 



Various ratios between the CMSSM and the SM of branching fractions of the main Higgs 
decay channels have been calculated with HDECAY 4.41 [16] for the regions prefered by the fit 
in the parameter space (see Fig. One notices an enhancement of the h ^ bb channel and a 
decrease of h — t+t^ by regards to the SM. Such a sensitivity makes potential measurements of 
the branching fractions attractive to discover a deviation from the SM and to determine the model 
parameters, even for SUSY mass scales beyond the LHC reach at y/s = 7 TeV or 8 TeV 



5. Conclusion 



We presented a global frequentist fit of the CMSSM and NUHMl parameter spaces, including 
up-to-date measurements in the flavor and electroweak sectors, the muon anomaly, astrophysical 
observations, the direct searches of supersymmetry at LHC and a Higgs mass of m/, = 126 GeV. 
The cuixent LHC exclusion leads to a low goodness-of-fit within the CMSSM, which worsens when 
requiring a Higgs mass above 125 GeV. The fit quality is increased in the non-minimal models 
NUHMl, despite a remaining tension due to the high correlations between observables in such 
model. The measurements of the Higgs branching fractions of the channels h — bb and h — t+t^ 
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Figure 6: Left: Parameter distributions for the LHC+m/, fit of the NUHMl with the 1 -dimensional lc7 in 
red and the 2-dimensional 2(7 in blue, and the best fit point marked by a star Right: predicted distribution 
of sparticle and Higgs boson masses from the LHC + w/, fit of the NUHM 1 . 



has the potential to indicate the deviation from the SM, even for SUSY mass scales beyond the 
LHC reach. 
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